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Abstract. The main objective of this article is to study the three-dimensional 
Rayleigh-Benard convection in a rectangular domain from a pattern formation 
perspective. It is well known that as the Rayleigh number crosses a criti- 
cal threshold, the system undergoes a Type-I transition, characterized by an 
attractor bifurcation. The bifurcated attractor is an (m — l)-dimensional ho- 
mological sphere where m is the multiplicity of the first critical eigenvalue. 
When m = 1, the structure of this attractor is trivial. When m = 2, it is 
known that the bifurcated attractor consists of steady states and their con- 
necting heteroclinic orbits. The main focus of this article is then on the pattern 
selection mechanism and stability of rolls, rectangles and mixed modes (includ- 
ing hexagons) for the case where m = 2. We derive in particular a complete 
classification of all transition scenarios, determining the patterns of the bifur- 
cated steady states, their stabilities and the basin of attraction of the stable 
ones. The theoretical results lead to interesting physical conclusions, which 
are in agreement with known experimental results. For example, it is shown 
in this article that only the pure modes are stable whereas the mixed modes 
are unstable. 



1. Introduction 

Over the years, the Rayleigh-Benard convection problem, together with the Tay- 
lor problem, has become one of the paradigms for studying nonequilibrium phase 
transitions and pattern formation in nonlinear sciences. There is an extensive lit- 
erature on the subject; see e.g. reviews by Busse [T], Chandrasekhar j^, Cross 
[3], Gethng [5], Koschmieder [6], Lappa [7], Ma and Wang [8], and the references 
therein. 

The problem is complete from the dynamic transition perspective (Ma and 
Wang [nun])- The main result in this direction is that the system always under- 
goes a Type-I (continuous) transition as the instability driving mechanism, namely 
Rayleigh number, crosses a critical threshold i?c, thanks to the symmetry of the 
linear operator, properties of the nonlinearity and asymptotic stability of the basic 
state at the critical threshold. Moreover, the system has a bifurcated attractor 
which is an (to — 1) -dimensional homological sphere where m is the number of 
critical eigenvalues of the linear operator. 

The main objective of this paper is to study pattern formation and the structure 
of the bifurcated attractor for the Rayleigh-Benard convection. The structure of 
the bifurcated attractor is trivial when m = 1. Namely, the attractor consists of 
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two attracting steady states approximated by the critical mode with opposite flow 
orientations. 

When m = 2, the picture is far from being complete. There are some known 
characteristics of this attractor such as the attractor must be homeomorphic to 
which contains either four or eight steady states connected by heteroclinic orbits 
or is a circle of steady states. 

From a pattern formation point of view, there is enough motivation to study the 
case m = 2. To understand the relative stabilities of steady states (patterns) to 
perturbations of other pattern types, there must exist at least two critical modes. 

In general the relation between the two horizontal length scales, for which m — 2, 
is nonlinear and hence it is very difficult to give a general characterization of every 
possible transition scenario. In this work, under the assumption that the wave num- 
bers of the critical modes are equal, we are able to give a complete characterization 
for m = 2 case. 

Depending on its horizontal wave indices ix and iy, a single critical mode can 
be either a roll (when at least one of ix or iy is zero) or a rectangle (when both ix 
and iy are non-zero) where ix and iy are non-negative integers which cannot vanish 
together. Thus there are three possible cases. Namely, (a) one of the critical modes 
is a roll while the other one is a rectangle, (b) both critical modes are rolls, (c) 
both critical modes are rectangles. 

In each case, we explicitly find nondimensional numbers which determine the 
number, patterns and the stabilities of the bifurcated steady states. We also deter- 
mine the basin of attraction of each of the stable steady states. 

In all the scenarios, we found that after the transition, only pure modes (rolls 
or rectangles) are stable and the mixed modes are unstable. Our result is conclu- 
sive when one of the critical modes is a roll type. When both critical modes are 
rectangles, we only have computational evidence. 

When both critical modes are rolls, the stable steady states after the transition 
are rolls. When both critical modes are rectangles, computational evidence suggests 
that the stable steady states after the transition are rectangles. When one critical 
mode is a roll and the other one is a rectangle, the stable states after the transition 
can be cither only rolls or both rolls and rectangles. 

The problem is usually studied in the infinitely extended horizontal domain 
setting which eliminates the effects of the boundaries in the horizontal directions. 
Our setting is a 3D rectangular domain with free-slip boundary conditions for the 
velocity, that is the fluid can not cross the boundaries but is allowed to slip. The 
thermal boundary conditions are adiabatically isolated side walls so that no heat is 
transferred through them and perfectly conducting top and bottom boundaries. 

Technically, the analysis is carried out using the dynamical transition theory (Ma 
and Wang [8j). One key ingredient is the reduction of the original system to the 
center manifold generated by the two unstable modes. The only modification that 
has been made is (following Sengul and Wang [H]), we expand the center manifold 
using a basis which differs from the eigenfunctions of the original linear operator. 
This allows us to passby the difficulties associated with determining the eigenpairs 
in terms of the system parameters. We also make use of computer assistance, 
namely a Mathematica code, which carries out numerous integrations which are 
due to the interactions of the critical modes with the non-critical ones. 
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The paper is organized as follows: In Section 2, the governing equations and 
the functional setting of the problem is introduced. Section 3 deals with the linear 
theory. We present our main results in Section 4. The proof of these theorems are 
given in Section 5. In section 6, we present the physical conclusions derived from 
our theorems. Finally, Section 7 is the conclusion section. 

2. Governing Equations and The Functional Setting 

With the Boussinesq approximation, the non-dimensional equations governing 
the motion and the states of the Rayleigh-Benard convection in a nondimcnsional 
rectangular domain O = {0,Li) x (0,^2) x (0, 1) C M'^ are given as follows; see [5] 
among others: 

+ (u • V) u = Pr (-Vp + Au + ROk) , 

09 

(2.1) — + {u-\7)e = w + Ae, 

Vu = 0, 

u(0) = uo, 9(0) = 00. 

The unknown functions are the velocity u = {u,v,w), the temperature 0, and 
the pressure p. These unknowns represent a deviation from a motionless state basic 
steady state with a constant positive vertical temperature gradient. In addition k 
stands for the unit vector in the z-direction. 



The non-dimensional numbers in (2.1) are the Rayleigh number R which is the 



control parameter and Pr, the Prandtl number. 

The above system is supplemented with a set of boundary conditions. We use 
the free-slip boundary conditions for the velocity on all the boundaries. Thermally, 
the top and the bottom boundaries are assumed to be perfectly conducting and the 
horizontal boundaries are adiabatically isolated. Namely, the boundary conditions 
are as follows: 

dv dw do 
u= — = — = — =0 atx = 0, Li, 
ox ox ox 

du dw do 

(2-2) dy='^-iry=dy=' ^^^^^'^^ 

du dv 

— ^ — ^0 ^0 atz = 0,l. 
dz dz 

For the functional setting, we define the relevant function spaces: 
H = {{u,0) e {n,R^) : \/ ■ u ^ 0,u ■ n \on^ O} , 
^^■^^ i/i ^ {{u,0) e H^n,R'') -.W ■u = 0,u-n\en=0,0 \,=o,i= O} . 

For = (u, 0),\et G : Hi^ H and Lr: H^^ H he defined by: 
LR(j) = (Pr7'(Au + ROk),w + AO), 
^^■^^ G(</.) = -(7'(u.V)u,(u.V)0), 

with V denoting the Leray projection onto the divergence-free vectors. The equa- 



tions J2.1| and (2.2 1 can be put into the following functional form: 



(2.5) -^=LRcj, + G{^), m^^o 
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The results concerning existence and uniqueness of ( |2.5[ ) are classical and we refer 
the interested readers to Foias, Manley and Temam [1] for details. In particular, 
we can define a semigroup: 

Sit) : 00 ^ 

Finally we define the following trilinear forms which will be used in the proof of 
the main theorems: 

G(0i, 02, 03) = - / (Ui • V)U2 • U3 - / (Ui • V)02 ' ^3, 

(2.6) Jn Jn 

Gs(01,02,03) = G(01,02,03) + G(02, 01,03)- 

3. Linear Theory 
We recall in this section the well-known linear theory of the problem. 

3.1. Linear Eigenvalue Problem. We first study the eigenvalue problem: 

Pr( Au + R9k-Vp) = (3u, 
(3.1) w + Ad^/3e, 

divu — 0, 



with the boundary conditions (2.2 1. Thanks to the boundary conditions, we can 
represent the solutions 0^ = (u5,0s), U5 = {us,vs,ws) by the separation of 
variables: 

us — Us sin(L^^s^7rx) cos(i^"'"Sj^7rj/) cos(sz7rz), 

vs = Vs cos{L'^^SxTrx)sm{L2^Sy7ry) cos{szTtz), 

(3.2) ^ ^ 

ws — WsCOs{L^ SjjTTx) cos(L^ Syiry) sm{szT^z), 

6s ~ 65 cos{L^^Sr^nx) cos{L2^ SyTiy) sin(s27rz), 

for S = {sx,Sy,Sz) where Sx > 0, Sy > 0, Sz > 0. It is easy to see that only 
eigenvalues l3s, S £ Z can become positive where 

Z = {{sx-, Sy, Sz) \ Sx>0,Sy> 0, Sz > 0, {sx, Sy) ^ (0, 0) and Sz ^ 0}. 

For S — {sx, Sy, Sz) G 2, the amplitudes of the horizontal velocity field can be 
found as: 



Us = — Y^s, Vs = — p Y^s- 



We define as, the horizontal wave number and 75, the full wave number by: 



(3.3) as = x^ + ^, IS 



sln^ / 52,^2 s27r2 



Huation 
Ap^ R 



Taking the divergence of the first equation in (3.1 1 we find: 

oz 



Now taking the Laplacian of the first equation, replacing Ap by the above relation 
and using (3.2) we obtain: 

7s(Pr7| + /3)Vl^5 - i?Pr40s = 0, 

(3.4) 

Ws-hl+ /3)0s = 0. 
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Figure 1. The selection of critical horizontal wave indices in the 
L1-L2 plane. The same coloring indicates equal wave indices. 



For each S G Z, the above equations have two solutions > /3| which satisfy the 
following equation: 

(3.5) 7l(7i + /5)(Pr7l + /3)-i?Pra| = 0. 

We find that amplitudes of the normalized critical eigenvectors as: 

(3.6) Ws^PUR)+^h = 



Now solving (3.5) for R at f3 — 0, the critical Rayleigh number can be defined as: 

(3.7) Rc = mini?s, i?s := 

^ ^ sez ' al 



From (3.7), one sees that for a minimizcr S — {s^, Sy, s^) of Rg, the vertical index 
S2 is 1. We will denote the set of critical indices S minimizing (|3.7|) by C: 



C = {S = (s,, sy, l)eZ\Rs<Rs', V5' e Z} 

For small length scale region, the map in Figure [l] shows the horizontal critical wave 
indices that are picked by the selection mechanism (3.7). 
It is well known that we have the following PES condition: 

r <0, \<Rc, 

(3.8) /3^(i?) = J = 0, A = i?„ yS eC, 

[>0, X>Rc, 

(3.9) Re/3(i?,)<0, V/3 ^ | 5 e C}. 



By (3.5 1, corresponding to S = {Sx,Sy,Sz) S Z, there are two eigenvalues I3g, and 
two corresponding eigenfunctions (p^g, i — 1,2. If a critical mode has wave index / 
then the corresponding eigenfunction is (j)] which we will simply denote by 0/. 
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3L2/2j, 




(a) A roll pattern with 
wave index J = {0,jy, 1). 

Figure 2. Regions of positive 
rectangular and a roll mode at 




Z.,/2i, 3Z.i/2i« £,, 

(b) A rectangle pattern 
with wave index / = 



iix,iy, 1). 

I negative vertical velocity of a 
mid-plane z = 1/2. 



Depending on the horizontal wave indices, there are two types of critical modes 
corresponding to two different patterns. If the wave index I — {ix,iy, 1) of a critical 
mode is such that one of the horizontal wave indices ix , iy is zero, the corresponding 
eigenfunction has a roll pattern. When both horizontal indices are non-zero, the 
corresponding eigenfunction has a rectangular pattern. Figure |2] shows a sketch of 
these patterns. 

3.2. Estimation of the critical wave number. As it will be shown, the dynamic 
transitions depend on the critical wavenumber a. In the case of infinite horizontal 
domains, the critical wave number is found to be a = 7r/-\/2 « 2.22 corresponding to 
a critical Rayleigh number Rc = 277r'*/4 w 658. For rectangular domains, the wave 
number is not constant and is a function of the length scales. Following estimates 
will be important in the physical remarks section: 



Lemma 1. Let a be the critical wave number minimizing (3.7 1. Then 



TT 



'^^ 21/3(22/3 + 1)1/2 - ^-55, forallLuL2, 

22/3^ 



a< , =g^3.10, i]Lx > 21/3^1 + 22/3^2.03, 

Vl -f 22/3 

a -> —=, Li — > oo. 

Proof. To estimate the dependence of the wave number a on the length scales Li, 
1/2, we define: 

L{m) ^ ((m + l)m)i/3((^ + i)2/3 ^ ^^2/3)1/2^ to e Z, m > 0. 

The sequence L{m) gives those length scales of Li for which the wave index changes 
assuming L2 is sufficiently small. As shown in Sengul and Wang |llj . when L{m — 
1) < Li < L[m) for some to > 1, we have the following bound on the critical wave 
number: 

TOTT TOTT 

(3.10) -— < a < 



L(to) L(to - 1) 

In particular, 

IT IT 

" - L(l) " 21/3(22/3 + 1)1/2^ 
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2.03 5 



Figure 3. The shaded region shows the possible values for the 
critical wave number a for a given Li. 



and 

27T 22/37r 
a < -rj— = , , if Li > L(l). 

- L{1) Vl + 22/3 ^ ' 



Finally noticing that, 



mir TT (m + l)7r tt 

— — — > —=, as m — >■ oo, 



L{m) 72' L{m) ^ 

we find that 



□ 



The bounds on the critical wave number as a function of the length scale Li 
which is obtained from (3.10) is shown in Figure [3j 

4. Dynamic Transitions and Pattern Selection 

We study the case where two eigenvalues with indices / = {ix,iy,l) and J = 
{jx,jy,l) are the first critical eigenvalues. This means that aj and aj minimize 
(3.7), thus the PES conditions (3.8), (3.9) are satisfied with C = {I, J}. The crucial 
assumption is that the corresponding wave numbers are equal, i.e. 

a ~ aj ^ OLj. 

Since I ^ J , without loss of generality, we can assume that > jx which ensures 
that iy < jy. By (3.3), we must have the following linear relation between the 
length scales: 

Thus two critical eigenmodes are possible only when Li and L2 lie on a line ema- 
nating from the origin in Figure [4] There are three possible cases depending on the 
structure of the critical eigenmodes, which are completely described by our main 
theorems. 

(a) a rectangle and a roll mode respectively (described by Theorem [T]), 

(b) both roll modes (described by Theorem [2]), 
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(c) both rectangle modes (described by Theorem [3|. 

These possible cases are illustrated by Figure |4] in the small length scale regime. 

Before presenting our results, we first summarize some of the known results which 
applies for the above setting; see Ma and Wang ^ [TU] : 

i) As the Rayleigh number R crosses Rc, the system undergoes a Type-I 
(continuous) transition. 

ii) There is an attractor E/j bifurcated on R > Rc such that for any (j>Q e H\r, 

dist{4)Q, Yir) — 0, as t ^ oo, 

where T is the stable manifold oi (p — Q with codim = 2. 

iii) Sfl is homeomorphic to and comprises steady states and the heteroclinic 
orbits connecting these steady states. 

iv) There are four or eight bifurcated steady states. Half of the bifurcated 
steady states are minimal attractors and the rest are saddle points. 

The dynamic transitions depend on the following positive parameter which in turn 
is a function of the parameters Pr, Li and L2'. 



(4.1) «5 = 



SPr^a^, 5 =(0,0, 2), 

^S#(i£+2Pr7^ + Pr2ii^), 5^(0,0,2). 



7r^(4a^- 
{Rs 
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i/rt (Rail) 



^8(Mixe^ ^ "fc (Mixed) 

*i(Rec)T _^ ^7 

• « • 

^1 (Mixedl ' ^6 (Mixed) 

Ifr3 (RoU) 



i/i, (RoU) 





•111 


l(Sec) 


*i'\(Rec) 


r y 





;i3 (Roll) 



(a) When a < c. 



(b) When c < a. 



Figure 5. The structure of the attractor after the transition R > 
Rc for different parameter regions when the first two critical modes 
are a roll and a rectangle. 



Here, a is the critical wave number, 7^ = + tt^ and Rc = jo? is the critical 
Rayleigh number. Moreover, for S = {sx, Sy, Sz), {sx,Sy) ^ (0,0): 



Li 



Also we let 



7^ 



Pra^ 



^ (Pr+1)74- 

4.1. One of the critical modes is a roll, the other is a rectangle. We first 
consider the case where an eigenmode with a roll structure and an eigenmode with 
a rectangle structure are the first critical eigenmodes. 

Theorem 1. Assume that I = (ix,iy,l) and J = (0,jj,,l), (ix > jy > iy > ^) 
are the first critical indices with identical wave numbers, aj — aj. Consider the 
following numbers: 

a = 'to, 0,2 + K2i^,0,2 + KQ,2iy,2, 

(4.2) b = Ko^o,2, 

For R > Rc, let us define the following steady state solutions: 



= g^/R - Re{X,cl)i + Y,<Pj) + o{{R - Rc)^'^), i = 1 

where 

X, = (-i)^a-i/2, = 0, 

= (-f)'(26)-l/^ 



X, =0, 



1,2, 



X, 



ab' 



ab' 



y^ = (-1)' 



Y^ = (-1)' 



2{c? - ab) ' 

c — a 
2(c2 - ab) ' 



(rectangle pattern) 
3, 4, ('roZ/ pattern) 

5,6, (mixed pattern) 

7,8, (mixed pattern) 



There are two possible transition scenarios: 
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i) If a < c then the topological structure of the system after the transition is 
Figure \5(a) ■ In particul 



as m 



ar: 

1) T,fi contains eight steady states ■(/'i, i = 1, . . . , 8. 

2) -01, ip2 (rectangles) and ^j^, ?/'4 (rolls) are minimal attractors of Tir, 

TpQ, ips (mixed) are unstable. 

3) There is a neighborhood U \ T of where T is the stable manifold of 
such that U — Uf^iUi with Ui pairwise disjoint and Ui is the basin of 
attraction of ipi, i = 1, . . . ,4,. 

4) The projection ofUi onto the space spanned by 4>i, 4>j is approximately 
a sectorial region given by: 

Ui n {X(j)i + Y(f>j I uj,s < arg{X, Y) < ^,,2}, i = 1, . . . , 4, 

Wi,i=7r-w, a;i,2 = 7r + a;, W24 = — 1^, ^2,2—^, 

Ws.i = n + uj, UJ3^2 = 27r — w, W4.1 — uj, W4.2 = tt — w, 

where 



(4.3) Lj = arctan , , , 

ii) If c < a then the topological structure of the system after the transition is 
as in Figure \5(b)\ In particular: 

a) T,fi contains four steady states ipi, i — 1, . . . ,A. 

b) ip3, V4 (rolls) are minimal attractors ofYiR whereas the ipi, 7/'2 (rect- 
angles) are unstable steady states. 

c) There is a neighborhood U \ T of where T is the stable manifold of 
such that U — uf^^tii with Ui pairwise disjoint and Ui is the basin of 
attraction ofipi, i = 3,4. 

d) The projection of Ui onto the space spanned by (pi, (j)j is a sectorial 
region given by: 

Ui n {X(l)i + Y(t>j I uj,^i < arg{X, Y) < ^,.2}, « = 3, 4, 

W34 = TT, UJ3^2 = 27r, UJi^l = 0, W4^2 = TT. 

Remark 1. In the special case jy = 2iy, the mixed solution corresponds to a regular 
hexagonal pattern. In this case we find a < c, hence the first scenario in Theorem^ 
is valid; see Remark^ 

4.2. The first two critical modes are both rolls. In this section we consider 
two critical modes both having a roll structure. Under the assumption that the 
wave numbers are equal, one of the rolls has to be aligned in the x-direction and 
the other one aligned in the y-direction. 

Theorem 2. Assume that I — {ix,Q,l) and J — (0,jj,,l) (ix > 1, Jj^ > are the 
first critical indices with identical wave numbers, aj = aj. Consider the following 
numbers: 

b — 2ko.o,2, 

(4.4) 

d = 2ko,o,2 + ^l^i^,jy,2- 

For R > Rc, we define: 



= g^R - R,{X,,Pi + Y,(Pj) + o{{R - R^)^'^), i = 1, 
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^ (Mixed) 



ipi (Roll) 

V 

^8 (Mix^'' 

iAi(Roll)T ^ ^ J 

f^ (Mixra), ^ ,^^6 (Mixed) 
^3 (Roll) 



Figure 6. The structure of the attractor after the transition R > 
Rc when the first two critical modes are both roll type. 



where 

X, = Y, = Q, i = 1, 2, (roll pattern) 

X, = 0, = (-l)^6"^/^ j = 3,4, (roll pattern) 

= (6 + c^)"^/^ = (-l)'Xi, i = 5, 6, (^mixed pattern) 

X, = -{b + dy^/^, Y, = {-lyX,, i = 7, 8, (mixed pattern) 

Then the topological structure of the system after the transition is as in Figure 
In particular: 

1) Sfl, contains eight steady states ipi, i ~ 1, . . . ,8. 

2) ipi, 'ip2, "037 "04 (rolls) are minimal attractors ofSn, ip^, ipQ, ipg (mixed) 
are unstable. 

3) There is a neighborhood U\T of where T is the stable manifold of such 
that U — Uf^iUi with Ui pairwise disjoint and Ui is the basin of attraction 
ofi/Ji, i = 1, . . . ,4. 

4) The projection of Ui onto the space spanned by (pi, <f>j is approximately a 
sectorial region given by: 

U, n {X(t>i + YPj I < arg{X, Y) < ^,,2}, i = I, . . . A, 

Wi,! = 37r/4, wi^2 = 57r/4, ^2,1 = -7r/4, W2,2 = 7r/4, 

^3,1 = 57r/4, a;3^2 = 77r/4, ^4^1 = 7r/4, ^4^2 = 37r/4. 

4.3. The first two critical modes are both rectangles. In this section we 
consider two critical modes both having a rectangular pattern with equal wave 
numbers, aj ^ aj. 

Theorem 3. Assume that I = {i^,iy,l), J = {jx,jy,l) (ix 0, iy 0; jx 7^ 0, 
jy ^ 0, ix 7^ jx, iy 7^ jy) are the first critical indices with identical wave numbers, 
aj — aj. Consider the following numbers: 
(4.5) _ 

a — ^0,0,2 + ^21^,0,2 + l^0,2iy,2, 
f — 'to, 0,2 + H2j^fl.2 + Ko,2j„,2- 
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. (Mixed) 
< 1 



^4 <Rec) 

(Ti (Rec) ^ ^ 

(Mixed) ^ (Mixed) 
^3 (Rec) 



1*4 (Rec) 

1*8 (Mixe^ (Mixed) 

((ri(Rec)i ^ ^ t 

(Mixe^,^^ ' (Mixed) 
(ia (Rec) 



(a) When a < e and / < e. 





/ i 

1 ^ 


>- < 


5*1 (Rec) T 





iln (Rec) 



(c) When f < e < a. 



(b) When e < a and e < /. 



1*1 (Rec) 




(friCRec) 



(d) When a<e< f. 



Figure 7. The structure of the attractor after the transition R > 
Rc for different parameter regions when the first two critical modes 
are both rectangle type. 



For R > Rc, let us define the following steady state solutions: 

V'z = gVR - Rc{X,(l)i + Y,(j}j) + o((i? - Rcf'^), z = 1, . . . , 8, 

where 

X, = {-lTa-^'\ = 0, 



X, = 0, 



X,, = 



X, 



I e-f ^ 
e2 - af ' 

[JET 



= i-l) 



1)7" 


-1/2 






/ e 


— a 


2(e2 


-afY 




/ e 


— a 


2(e2 


-afY 



i = f , 2, (rectangle pattern) 

z = 3, 4, (rectangle pattern) 

z = 5, 6, (mixed pattern) 

i = 7,8, (mixed pattern) 



There are four possible transition scenarios: 

i) If a < e and / < e then the topological structure of the system after the 
transition is as in Figure \7(a) . In particular: 
I) T,fi contains eight steady states ipi, i — 1, ... ,8. 
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2) ipi, -02, ip3, ip4 (rectangles) are minimal attractors ofJ^n, ip^, ipQ, ipj, 
ipg (mixed) are unstable. 

3) There is a neighborhood U \ T of where T is the stable manifold of 
such that U — uf^iUi with Ui pairwise disjoint and lAi is the basin of 
attraction of tpi, i — 1, . . . , 4. 

4) The projection oflAi onto the space spanned by (pi, <f>j is approximately 
a sectorial region given by: 

Ui n {X(t>i + Y(t)j I < arg{X, Y) < uj.^}, i = 1, . . . , 4, 

Wia=7r-w, a;i,2 = vr + a;, ^24^-0;, ^2^2=^^, 

0-13.1 = TT + cij, a;3.2 = 27r — w, ^4.1 = uj, 0^4.2 — t: — uj, 

where uj — arctan , 

ii) If e < a and e < f then the topological structure of the system after the 
transition is as in Figure 1 7(b )\ In particular: 

a) Y,f{ contains eight steady states ^Ji, i = 1, . . . , 8. 

b) '05, TpQ, ipf, 08 (mixed) are minimal attractors ofUji whereas ipi, "02, 
"03, -04 (rectangles) are unstable. 

c) There is a neighborhood IA\T of where T is the stable manifold of 
such that U — Uf^^Ui with Ui pairwise disjoint and Ui is the basin of 
attraction 0/ -01, * = 5, . . . , 8. 

d) The projection of Ui onto the space spanned by 4>j, 0j is a sectorial 
region given by: 

U, n {X(j)i + Y4,j I < arg{X, Y) < ^,,2}, i = 5, . . . , 8, 

W5,i = 0, a;5,2 = 7r/2, wg.i = Sn/2, ^5,2 = Stt, 

W7,i = TT, a;7,2 = 37r/2, cjg.i = 7r/2, = tt. 

iii) If f < e < a then the topological structure of the system after the transition 



is as in Figure 7(c) In particular: 

a) Sfl contains four steady states ^i, i = 1, . . . , 4. 

b) "011 "^2 (rectangles) are minimal attractors ofYiR whereas 03, "04 (rect- 
angles) are unstable steady states. 

c) There is a neighborhood U\T of where T is the stable manifold of 
such that U — uf^iUi with Ui pairwise disjoint and Ui is the basin of 
attraction ofipi, i = 1,2. 

d) The projection of Ui onto the space spanned by 0/, 0j is a sectorial 
region given by: 

U, n {X0/ + r0j I < arg{X, Y) < ^,,2}, « = 1, 2, 

wi.i = 7r/2, 1LJ1.2 = 37r/2, ^2,1 = -7^/2, ^2,2 = 7r/2. 

iv) If a < e < f then the topological structure of the system after the transition 
is as in Figure \7(d)\ In particular: 

a) Sjj contains four steady states ^i, i — 1, ... ,4. 

b) Tp3, "04 (rectangles) are minimal attractors ofSji whereas the ipi, ^2 
(rectangles) are unstable steady states. 
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c) There is a neighborhood U \ T of where T is the stable manifold of 
such that U — Uf^^Ui with lAi pairwise disjoint and Ui is the basin of 
attraction oftpi, i = 3,4. 

d) The projection of Ui onto the space spanned by is a sectorial 
region given by: 

U, n {X(j)i + Y(j)j I < arg{X, Y) < ^^,,2}, « = 3, 4, 

= TT, UJ3^2 = 27r, tJ4,l = 0, W4^2 = TT. 

5. Proof of the Main Theorems 

First we give the preliminary setting that wih be used in the proof of the main 
theorems. 

The first step is to find the adjoint critical eigenvectors. The adjoint equation of 
(fsll) is: 



(5.1) 



Pr(Au* - Vp* 

A0* +RPtw* 



The eigenfunctions of (5.1) can be represented by the separation of variables (3.2). 
Also the eigenvalues of (5.1) are same as the eigenvalues of (3.1), i.e. /3 satisfies 

Ijoi: 

e*s 



(3.5). We find the amplitudes of the critical adjoint eigenvectors as: 
(5.2) W*s=f3l{R)+^l 



RPt, 



where PgiR) satisfies the PES condition (3.8). 

Let / and J be the indices of the critical modes, i.e. C — {I, J} in (3. 
denote: 

0, = 0,7 = f5{R) = f3]{R)^^]iR). 
We study the dynamics on the center manifold, we write: 

= 0= + 



We will 



in (2.5), where $ is the center manifold function and 

0= = xi0/ + a;20j. 



Multiplying the governing evolution equation (2.5) by the adjoint eigenvectors 
0} and we see that the amplitudes of the critical modes satisfy the following 
equations 

dxi 1 

(5.3) 

(G(0),0}) 



1 



(G(0),01) 



The pairing (•, •) denotes the L^-inner product over il. 



By ([3^ and ( |3.5| ), 
(5.4) P{R) ^ g{R- Rc) + oiiR- Rc)^), as R 

with 

Pra^ 



^=(p7TW^"- 

We write the phase space as 

H = Ei®E2, -El = span{07,0j}, 



E2 — E^. 
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When the hnear part of (5.3) is diagonal, we have the fohowing approximation 
of the center manifold; see Ma and Wang [5]: 

(5.5) = P2G(0=)+o(2), 

where Cr ~ Lji j^^, P2 is the projection onto E2 and 

o(n) = o(N") + 0(N"|/^(i?)|). 



Using (5.5) we can write the center manifold function as: 
y) = a;^$i + xy<^2 + y^*3 + o(2) 



Using this approximation and the bilinearity of G, we can write (5.3) as: 
dx 

dy 



(5.6) 



: l3{R)x + x{aix'^ + a2y^) + 
l3iR)y + y{bix^ + b2y^) + 



o(3), 
0(3). 



In dynamic transition problems, the center manifold is generally expanded using 
the eigenfunctions of the original linear operator. However, following Sengul and 
Wang [TT], we will expand the center manifold using a different basis. Namely we 
will consider the eigenfunctions U5 of the Stokes equation for the velocity together 
with eigenfunctions 9s of the Laplace equation. The main advantage of such an 
expansion is that the eigenvalues and eigenfunctions are independent of the system 
parameters, namely the Prandtl number Pr and the Rayleigh number R, while still 
spanning the same functional space (2.3) which leads to computational advantages. 

For this reason we turn to the following eigenvalue problem with the boundary 
conditions (2.2) of the original problem: 

Aus - Vp ^ pus, 
A9s = pes, 
divug = 0. 

By the classical theory of elliptic operators, the eigenvectors {e^ = (u5,0), = 
(0, 0s)} form a basis of the phase space H. Moreover, eg can be expressed by the 
same separation of variables (3.2 1. There are three cases to be considered. 
• If {sx, Sy) — (0,0) and 7^ 0, then = and 

{0,0s), es = i. 



el 



If, 



Sy^O and = 0, then there are eigenmodes which have the form 



(u, 0) with u 



{u, V, 0). For such modes, it can be easily verified that: 
{Gir),e)=0. 



Thus by (5.51, such modes will not be present in the lowest order approxi- 
mation of the center manifold function. 

If + 7^ and 7^ 0, then the multiplicity of an eigenvalue is two and 
the eigenvectors are: 



e's 



(U5,0), 
(0,^5), 



Ws 
Qs-- 



1. 



The following lemma guarantees that the projection P2 in (5.5) can be ignored 
from a computational point of view. 
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Lemma 2. For i = 1,2, 

P2e's = e'c for S^C, 

(5.7) ^ s s ^ , 

{P2Gm,e's)^0 for See. 

Proof. First note that: 

El = span{4>}, cjjj} C span{e5, e| | 5 G C}. 

Thus 

E2=El D span{e^, e| | 5 ^ C}. 

Thus we have the first equation in (5.7). 
Since 

span{4, (j>l\SeC}^ span{e^, \ S & C}, 
there must exist constants cs^i ^ 0, 03,2 7^ such that 

^265 = C5a0|, P2e| = cs,2(/'|, for S eC. 
Since by direct computation we have 

(F2G((/)^),4) -0. for^eC, 



we also have second equation in (5.7). □ 
Now we write: 

(5.8) <^'si^,y)es + o{2), 

ses, 1=1,2 

where S denotes some index set which will be specified later. Here <i>g are quadratic 
polynomials in x and y. 
Let 

Z™" ^{K = (fc, 0, 1) or X (0, fc, 1) : fc 7^ 0, aK = a}, 
^rec ^{K={ki,k2,l):ki^Q,k2^0,aK^a}. 

Lemma 3. For Si = (0, 0, 2), we have: 
where for K e {/, J}, 



(5.9) $ 



K _ ) -^[0,1]"^, ^fKeZl^ 



Si 



Also for Ki e K2 € Z™"; 

(5.10) G,(0x,,es,>;,J = iG,(0K.,e5,,01^J = ^^^'^^'^^ [0, Pr]^- 

Proof. If (sa;, Sy) = (0, 0) and ^ 0, then by (5.5), for K e {/, J}, we have 
^5 j^-^^ ^ (G('/'g,0g),e|) ^ (G('/'j<-,0K),e|) 

Note that 

(5.12) (e|,e|)= / sin2 27rz=^^. 
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For Ki e Z^"^"^, K2 £ -Z™ , a direct computation yields: 



(5.13) 



[0,1] 



Here Gg is the trilinear operator defined in (2.6|. Now, (5.9) follows from (5.11), 
(5.12) and (5.13). Also, (5.10) follows from an easy computation. □ 



Lemma 4. If S ^ C, S ^ (s^, Sy, Sz), {sx, Sy) ^ (0, 0) and Sz 7^ 0, the 



(5.14) 

Here 

(5.15) 

and 

(5.16) 



$0 = 



$1 



-4, 



(G(</'^),4; 

(G(</.^),4; 



-PrRsJs RPr 
1 -7| 




(5.17) 
Proof. Let 
(5.18) 
That is: 

Let 



^5 



^6 



m, n = 1, 2. 



il 



2 '^x'^'^l 2 '^^/^■^2 9 
cos ; COS COS SzTTX^. 



Li 



Clearly is equal to the definition in (5.16). Also it is easy to compute that As 



is the same as given in (5.15). Now the lemma can be proved by taking the inner 
product of (5.5 ) by eg, i = 1, 2 and using Lemma [2] □ 



Notice that 



dety^s = vlPr{Rs-R). 



Since Rs > Rc for S ^ C, the determinant of A is always positive when R is close 
to Rc- This guarantees that (5.14) can be solved for $g and $|. 
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Finally note that for S = {sx, Sy, Sz) with {sx-, Sy) ^ (0, 0) and Sz ^ Q sX R = Rc 
we have: 



s,(l>*s) = / usu*g + vsvs + wsw*g + ^ 



(5.19) 



^M^{il + m)f + RP^) \r=r. 



===05(Pr + l)75- 
Notation 1. We will use the following notation: 

G[(pi,(pj,es) = 

G(0/>J,4) 

We also define the following indices which will be used through out the proofs. 
5i = (0,0,2), ^2 = (2i.,0,2), ^3 = i0,2iy,2), = (z,,i^ + j„2), 

Sb = {ixdv -*!;;2), 56 = {ixjy,2), Sj = {i^ + jx.iy +jy,2), 
S% = {ix - jx,iy + jy,2), Sg = {ix + jx-.-iy + jy,2), 

Sio - {tx-Jx,-iy+Jy,2), Sn = (2j,,0,2), Si2 = (0,2jj„2). 

5.1. Proof of Theorem [l| Now we assume that / = {ix,iy, 1) and J = iO,jy, 1), 
{ix > 1, jy > *a ^ 1) are the first critical indices with identical wave numbers, 
aj = aj. Using (5.5), we find that the lowest order approximation of the center 
manifold is spanned by the eigenvectors having indices 5*^ as given in (5.20): 

$(x,2/) = a;2($^^esi + ^ <^>ses) + xy ^ «'ses + y^^SiCSi ■ 

S— 82,33 S— 34,85 



(5.20) 



Here 



where according to (5.9) we have: 



The coefiicients in (5.6) are as follows: 



(5.21) 



ai 



a2 



hi 



|^$^^G,(07,e5,,01)+ ^ $5.G.(0z,e5.,01)j , 
-1— U^^G,(07,es,,01)+ 5] $5.G.(</.j,es,,0D , 



1 



i=4,5 



■$^^G,(0j,es,,<^}). 



All the terms involving Si above can be computed using the Lemma [3j 
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By direct computation: 
(5.22) Gs{4>i,es,(j)*j) ^ Gs{4>j,es,4>*i), 8 = 84,35. 



Putting ([5^, ( |5.22[ ) and ( |5.19| into ( |5.21[ ) we find that: 

(5.23) 02 = 2bi. 

By direct computation we can obtain the following: 



(5.24) 



G 



s\<Pli<P.J, 



(5.25) 

where 
(5.26) 



1 




a-2 












) = -Vs 








. 


1 








F 


'i-RcY 


1 








Pri?c7 



8 — 82, 83, 
8 = 5*4, S'5, 

, 8 = 82, 83, 
, 8 = 84, 85, 



Vs = 



32a2 



Using (5.11), (5.141 and (5.24) , we find: 



(5.27) 



$1 
$1 



$1 



$1 



1 7r(4a2 - a|)72 
16 Pr a'^{R^ - Rs) 

1 t:{W - Q!|)72 



Pr i?cQ;2 + 747I 
Pr + 1^ 



S* — 5*2, 5*3, 



5 — 54, 5*5 



4 Pr a4(i?,-i?5) 

Now putting (5.251, (5.27), ( 5.13[ ), (5.9) into (5.21) and normalizing the results by 
(5 28) 

^ ■ ' 2ioPra2^0/,0}) 27a2pi.(i + Pr)' 

and using (5.23), the equations ( |5.6[ ) become: 
dx 
~dt 



l3{R)x - x(ax^ + 2cy2) + o(3). 



(5.29) 



dt 



l3iR)y-yicx^ +2by^) + o{3), 



where a, b and c are as defined in (4.2). Note that Kg. > in (4.1) since R^ <^ Rsi ? 
i = 1, . . . , 5 and we have the following relations: 

(5.30) 0<6<a, 0<6<c. 



Now we consider the approximate steady state equations of (5.29): 

(3{R)x - x{ax^ + 2cy2) = 0, 
PiR)y - yicx^ + 2by^) = 0. 



(5.31) 
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Let us define 
(5.32) 

There are two cases to consider. 



2{c-b)' 



i) If c < a, tlicn tlie equations (5.31 1 have only four straight hue orbits on the 
lines y — and x — 0. And the following solutions of (5.31 ) are bifurcated 
on /3 > 0: 



(5.33) 



^4 



Y+ 



(0,±1 



ii) If c > a then there are four additional straight line orbits on the lines 
y = ztutx. Note that in this case by ( 5.30 1, — a6 > and there are four 



(5.34) 



additional solutions bifurcated on /3 > which are given by: 

^ =1,2. 



(-l)"(l.±")l//3 



Now the Jacobian of the vector field in (5.31) is: 



(5.35) 



J = 



(3 — 3ax^ — 2cy^ —4cyx 

—2cyx P — cx^ — 6by^ 

The Jordan canonical forms of the Jacobian matrix evaluated at the steady states 



in (5.33) and (5.34) are as follows 



(5.36) 




-1 



, j(r±) = 2/3 



fe-c 
26 











(b— c)(c— a) 
ab—c"^ 



The stability of the steady states can be found by using (5.36). Since c > h, Y±. 
are always stable on R > Re- If c > a then X± are stable and Z± are unstable 
on i? > Re- On the other hand if c < a then X± are unstable and Z± are not 
bifurcated on R > Re- Thus we have two transition scenarios and the results are 
shown in Figure [5j 



Remark 2. In the particular case jy — 2iy, the equations in (5.291 can he reduced 
further. In this case we have: ag^ = as^ and agg = as^ which implies that 
7^2 = ISi and 75, — 7s5 . This in turn implies — k-s^ and ks, — ks^ and we 
get the relation 2a — b + c. This implies c > a since c — a — a — b>0. 

5.2. Proof of Theorem [2j We point out the differences from the previous proof. 
We have / — {ix,0, 1) and J — {0,jy, 1) {ix > 1, jy > 1) as the first critical indices 
with identical wave numbers, aj = aj. First the center manifold is given by: 



^ix,y) = x^'^?'s,<Psi+xy<^SeeSe 



yHi^es,+oi2). 
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Using this approximation, the coefficients in (5.6) are as fohows: 

1 ^ , . 

ai 



(5.37) 



02 

bi 
b2 









b}) 


1 








{<pl,(, 




1 












1 





(<I>^^G3 65,, 0J) + $SeG.(0j, es,, ^J)) , 
<i>^,G,((/.,,, 65,, 



The coefficients ai and 62 are equal to 62 in the proof of the first theorem. Thus 
we only need to find 02 and bi. A quick computation shows that: 

(5.38) G,(</)/, 656>}) = G,(0j, 65„ 0J), 
which shows that: 

(5.39) 02 = &i. 

By direct computation we can obtain the following: 

-2 

(5.40) G,(0,,0j,65) = -2??5 " , S = S^ 



(5.41) 



Gs(0j,e5, (t)*j) = 775 



5* — Sf. 



7 



a"" 
Pr i?c7^'^ 

where rj is given by (5.26). Using ( |5.11[ ), (5.14) and (5.40) , we ffiid: 

Pr R^a^ + 7^7! 



(5.42) $5 = 



1>5 
$1 



1 7r(4Q2 _ Q,2 

2Pra4(i?,-i?5) 



S — Sp, 



By normalizing all terms in ( 5.37[ ) by the normalizing factor ( 5.28 ) and using ( 5.39 ), 
the equations (5.6 1 become: 

dx 

~di 
dy 
dt 



(5.43) 



(3{R)x - x{bx^ + dy^) + o(3), 
P{R)y-yidx^ + by^) + o{3). 



Here & > and d > are deffiiition (4.4). Since /t > 0, we always have b < d. Now 
we consider the approximate steady state equations of (5.43): 

l3{R)x - x{bx'^ + dy^) =0, 

PiR)y-y{dx^ + by^)^0. 



(5.44) 



The equations (5.44) have always eight straight line orbits on the lines y ~ 0, x — 
and y = ±x. The eight solutions of (5.44) which are bifurcated on /3 > are: 
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The Jordan canonical forms of the Jacobian matrix evaluated at the steady states 



in (5.45) are as follows: 

J{X±) = 2/3 



(5.46) 



J{Zl) ~ 2/3 



-1 

-1 







d+b 



, J{Y±) = 2(3 



, 1 = 1,2. 



b-d 
2b 





Using (5.461 we can find that, since b < d, X± and Y± are stable and Zj., i = 1,2 
are unstable on R > Rc . That finishes the proof. 

5.3. Proof of Theorem [3} Again we point out the differences from the previous 
proofs. We have / = {ix,iy,l), J = {j^,jy,l) {i^ ^ 0, iy ^ 0, ^ 0, jy ^ 0) 
as the first critical indices with identical wave numbers, a/ = aj. First the center 
manifold is given by: 

=a;2($^^(/)Si + ^ ^ses) + xy ^ $ses 

S— 82,33 S— Sr .Ss,Sq,S 10 

+ y\'^'k<Ps^+ *ses) + o(2). 

S=S2..S3 

Using this approximation, the coefficients in (|5.6|) are as follows: 

{(Tit. m"^) 

1=2,3 



(5.47) 



02 



bi 



1 








61) 


1 








('/'/,'; 


61) 


1 












1 













1=7,.. .,10 



i=7,...ao 



1=11,12 



Oi and 62 is computed in the same way as ai in the proof of the first theorem. So 
we will only find 02 and bi. Also using 

(5.48) Gs((/)/, eg, (/)}) = Gs(</)j,es,(/)J), S = Sj, Sg, Sg, Sio, 
we find that: 

(5.49) 02 = bi. 

By direct computation we can obtain the following: 

0,-2 



(5.50) 



Gs{(l)i,(j)j,es) = --VS 



1 



(5.51) Gs{<i}j,es,^*i) = -ris 



S — St, Ss, Sg, Sio, 



S — Sj, 5*8, Sg, Si 
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where ij is given by (5.26). Using (5.111, (5.14) and (5.50) , we find: 
(5.52) 



$1 



1 7r(4a2 _ Q,2 



Pr + 7^71 
Pr + 7^ 



S — Si, iSg, 5*9, S'lo- 



By normalizing all terms in ( 5.47[ ) by the normalizing factor ( 5.28 ) and using ( 5.49 ), 
the equations (5.6 1 become: 



(5.53) 



dx 
'dt 
dy 
dt 



: j3{R)x - x{ax'^ + ey"^) + o(3), 
:/3(i?)y-y(ea;2 + /y2)+o(3). 



Here a, e and / are positive numbers defined in (4.5). 



(5.54) 



Now we consider the approximate steady state equations of (5.53) 

l3{R)x - x{ax'^ + ey^) = 0, 
/3(i?)2/-y(ex2 + /y2) =0. 



The analysis of the equations (5.54) is similar to the analysis of the equations (5.31 ) 
given in the proof of Theorem [1] Thus we omit the details. 

6. Physical Remarks 

In this section we will use the main theorems to derive some physical conclusions. 

Before going into details, we have to make a remark about the critical wave 
number a since it is one of the parameters determining the transition numbers. 
Although a depends on the length scales, thanks to Lemma [T] we have certain 
bounds on its range of values. To recall, these are: 

i) a > 1.55 regardless of the length scale. 

ii) If one of the length scales is greater than 2.03 then a < 3.10. 

Taking a look into Figure one sees that the only case which is not covered by 
taking a in the range 1.55 < a < 3.10 is that of two critical rolls with indices 
/ = (1, 0, 1) and J — (0, 1, 1) which happens when Li = L2 < 1.69. 

There are only three possible cases when two modes with equal wave numbers 
become unstable simultaneously. Namely these two critical modes can be, a rect- 
angular and a roll mode, both roll modes, both rectangular modes. We investigate 
each case separately. 

6.1. The first two critical modes are a roll and a rectangle. We first consider 
two critical wave indices / = {ix, iy, 1) of a rectangular pattern and J = (0, jj,, 1) 
of a roll pattern with equal wave numbers a ^ aj = aj; see Figure^ We define 
the following number: 

ly 

Since we are assuming that aj = aj, we have 1 < iy < jy and A > 1. Notice that 
A is the number of rolls to the number of the rectangle columns in the direction of 
rolls. Figure shows that the possible values of A in the small length scale regime 
max{ii,L2} < 5.5 are 4/3, 3/2, 2, 3, 4. 

Now to use Theorem [1] to describe the pattern selection after the transition, we 
have to compute the transition numbers a and c given by (4.2). These numbers 
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depend on three parameters Pr, Li and L2. Equivalently one can use Pr, A and a 
as the parameters determining a, b and c. 
Using 



ct I (y. I (y, 

as^^2\la^ - —, as3=2-, as^ ^ J 2a{a + as, ^ J 2a{a ~ -), 



in (4.2 1, we can compute the numerical values of a and c for a given value of A, a 
and Pr. 

There are two cases that can happen, depending on whether a < c oi c < a. 
In i^z(7wre[Pl the regions where c < a is shown for several parameter regimes. The 



results show that both transition scenarios described by Theorem 2.1 are possible. 
In particular there are two parameter regimes such that c < a, namely when Pr < 
0.2 and 1 < A < 1.4 and for large A or large Pr. 

Now some remarks about the basin of attraction of the rolls and rectangles are 
in order. When a < c, both rectangles and rolls are stable but their basins of 
attraction which are sectorial regions depend on an angle ui which depends on Pr, 
a and A. In the particular case A = 2, the mixed modes have a hexagonal pattern 
and we find that a < c. Moreover, ui is independent of a and Pr and is found to be 
w = arctan 1/2 « 26.57°. So in this case the basin of attraction of rolls consists of 
two sectors each of which has an angle of tt — 2w « 126.87° while the basin for the 
rectangles have an angle of 2uj « 53.13°. This means that rolls will attract a wider 
region of initial conditions than rectangles do. 

6.2. The first two critical modes are both rolls. Now we consider two roll 
type critical modes. By the assumption of equal wave numbers, the rolls has to 
be perpendicular to each other, i.e. / = {ix,0,i) and J = (0,jj,,l). In this case 
after the first dynamic transition, always the rolls are stable and the mixed states 
are unstable. Moreover, rolls with index / and rolls with index J have uniform 
attraction basins. 

6.3. The first two critical modes are both rectangles. Now we consider the 
case where the first two critical modes both have rectangle patterns, i.e. / = 
{ix,iy,^) and J = {jx,iy, 1) {ix > jx > 1; ^ h — 1) the first critical wave 
indices. In this case, the dynamic transitions depend on the numbers a, e and / 
given by Theorem [3] To calculate these numbers we define the parameters A and 
B: 

A=^>1, B=^>1. 

"iy Jx 

Using the definition, we find that: 



2 A V ' „2 „2 4 ^ „2 

2«x.o B^A^ — 1 " B^A^ — 1 



2 A^{B'-1) 2 2 . ^ -1 2 

^2jx.O = 4 ^2^2 _ I " ' "0,2j„ = ^ g2^2,l " ' 



Cti ,1 / T\n^ ,A =2 [a H an,- n+A ttn 9i , 771,71=1,2. 

*x + (-1)'"Ji:(-1)"*!,+Jh y 4^ 2lx,() ^ 0,2ly ) ' ' ' 

Using this, we computed a, e and / for several choices of Pr, a, B and C. Our 
numerical calculations revealed that a < e and / < e for a vast amount of parameter 



choices. This means that the transition scenario is described as Figure 7(a) Then 



the rectangles with index I and J are both stable and the mixed modes are unstable. 
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However, the other transition scenarios can also be possible as we observed that 
f < e < a when the Prandtl number is small, one of ^ or is less than 2 and 
A ^ B. For an example see Figure In this case the transition scenario is 
described as Figure 7(c)\ Hence only rectangles with index / are stable and the 
rectangles with index J are unstable. 

In particular, our numerical investigations suggest that one or both of the pure 
modes (rectangles) are stable while the mixed modes are unstable. 



7. Conclusions 

In this paper, we discuss the dynamic transitions of Rayleigh Benard (RB) con- 
vection from a perspective of pattern formation. We focus on the case when two 
eigenvalues cross the imaginary axis simultaneously. This allows us to compare the 
stability of a pattern with respect to perturbations of other pattern types. Our 
main assumption is that the wave numbers of the critical modes are equal. Under 
this assumption, we classify all the possible transition scenarios and determine in 
each case the preferred patterns and their basins of attraction depending on the 
system parameters. 

The pattern of a simple critical mode is either a rectangle or a roll. Thus there are 
three possible cases when there are two critical modes: (a) one mode is rectangular, 
the other mode is a roll, (b) both modes are rolls, (c) both modes are rectangles. 

The following are some general characteristics of the transition for the RB con- 
vection with which already known (Ma and Wang [9l HO]): 

1) The transition is Type-I. In particular, there is an attractor bifurcating 
on i? > Re- 

2) Ylji is homeomorphic to 5*^ which comprises steady states and the connect- 
ing heteroclinic orbits. 

The following are the results due to our main theorems: 

3) In all the scenarios, we found that only pure modes (rolls or rectangles) 
are stable and the mixed modes are unstable. Our result is conclusive 
(analytical proof) when one of the critical modes is a roll type. When both 
critical modes are rectangles, we only have computational evidence. 

4) When both critical modes are rolls, the stable steady states after the tran- 
sition are rolls. When both critical modes are rectangles, computational 
evidence suggests that the stable steady states after the transition are rect- 
angles. When one critical mode is a roll and the other one is a rectangle, 
the stable states after the transition can be either only rolls or both rolls 
and rectangles. 

5) When both rolls and rectangles are stable after the transition, these states 
have non-uniform sectorial basin of attractions. In the particular case, 
where the mixed states have a regular hexagonal pattern, the angle of the 
sector for rolls is 126.87° while the angle of the sector for the rectangles is 
53.13°. Thus rolls attract a wider range of initial conditions, making them 
a more preferable type of pattern. 
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Figure 8. The selection of the horizontal wave indices (ixiiy)- 
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Figure 9. The shaded regions show the parameter regimes where 
c < a when two critical indices are / = {ix, iy, 1) and J = (0, jj,, 1). 
Here A = jy/iy 
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Figure 10. The parameters a — e and / — e in the small Prandtl 
number regime. 
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